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ABSTRACT 

Aims. Our aim is to investigate the history of mass assembly for galaxies of different stellar masses and types. 

Methods. We selected a mass-limited sample of 4048 objects from the VIMOS VLT Deep Survey (VVDS) in the redshift interval 0.5 < z < 1.3. 
We then used an empirical criterion, based on the amplitude of the 4000 A Balmer break (Dn4000), to separate the galaxy population into 
spectroscopically early- and late-type systems. The equivalent width of the [OII]3727 line is used as proxy for the star formation activity. We also 
derived a type-dependent stellar mass function in three redshift bins. 

Results. We discuss to what extent stellar mass drives galaxy evolution, showing for the first time the interplay between stellar ages and stellar 
masses over the past 8 Gyr. Low-mass galaxies have small Dn4000 and at increasing stellar mass, the galaxy distribution moves to higher Dn4000 
values as observed in the local Universe. As cosmic time goes by, we witness an increasing abundance of massive spectroscopically early-type 
systems at the expense of the late- type systems. This spectral transformation of late-type systems into old massive galaxies at lower redshift is 
[ a process started at early epochs (z > 1.3) and continuing efficiently down to the local Universe. This is also confirmed by the evolution of our 
. type-dependent stellar mass function. The underlying stellar ages of late-type galaxies apparently do not show evolution, most likely as a result of 

■ a continuous and efficient formation of new stars. All star formation activity indicators consistently point towards a star formation history peaked 

■ in the past for massive galaxies, with little or no residual star formation taking place in the most recent epochs. In contrast, most of the low-mass 
systems show just the opposite characteristics, with significant star formation present at all epochs. The activity and efficiency of forming stars 
are mechanisms that depend on galaxy stellar mass, and the stellar mass assembly becomes progressively less efficient in massive systems as time 
elapses. The concepts of star formation downsizing and mass assembly downsizing describe a single scenario that has a top-down evolutionary 

' pattern in how the star formation is quenched, as well as how the stellar mass is grown. The role of (dry) merging events seems to be only marginal 
' at z < 1.3, as our estimated efficiency in stellar mass assembly can possibly account for the progressive accumulation of observed passively 
evolving galaxies. 
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1. Introduction 

How galaxies form and evolve with cosmic time is one of the 
key questions in observational cosmology. A fundamental role 
towards answering this question is played by deep surveys sam- 
pling over a thousand galaxies on large portions of the sky. The 
global star formation history is reasonably well known out to 
very high redshifts an d with considerable details up to z ~ 1 (see 
for a recent summary ^ Hopkins & Beacomll2006l) . but the role of 
the various physical mechanisms contributing to the assembly 
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of galaxy stellar mass is still unclear, as is their importance at 
different epochs. 

In agreement with the first formulation introduced by 
ICowie et all (Il996l) . a downsizing scenario for galaxy for- 
mation persistently emerges in several observational studies 
(Brinchmann& Ellis 2000; Gav azzi et alj 120021; iKodama et al 
2004; Bauer etal. 2005; Feulner et al.' '20053'!?; 'Juneau et al 
2005; Bundy et al. 2006; Borch et al. 2006; Cimatti et al. 200^^ 
iCucciati et all 120061; iPozzetti et all 120071) . However, observed 
galactic-scale properties like the quenching of star formation 
activity in massive galaxies not yet h ave been reproduced 
well within the ACDM framework (e.g. iDe Lucia et al. 1 120061; 
Bower et al. 2006) that successfully describes the hierarchical 
growth of dark matter halos and the galaxy clustering properties. 

The observed discrepancies between models and observa- 
tions are probably due to our limited ability to reproduce the ac- 
tual physics on galactic scales with simple recipes. Perhaps some 
feedback mechanisms are still missing in models. Thus, a phe- 
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nomenological approach to the problem of the stellar mass as- 
sembly as a function of cosmic time, i.e. one that relies as much 
as possible on observational data, can suggest new implementa- 
tions in current semi-analytical models. 

In this paper we use the first epoch obse rvations of the 
VIMOS-VLT Deep Survey (VVDS. lLe Fevre e t al. 2005) to de- 
rive stellar masses, stellar ages and star formation efficiency for a 
statistically significant sample of galaxies. We analyse how stel- 
lar mass is assembled in active and quiescent galaxies classified 
following a spectroscopic scheme, and how this process evolves 
with cosmic time using spectroscopic diagnostics. This analysis 
provides constraints on the hierarchical scenario thought to gov- 
ern dark matter assembly and can be pursued using VVDS data 
starting from a time when the Universe had ~ 30% of its present 
age. 

The present work is organised as follows: the data and sam- 
ple selection are presented in Section 2; the methodology we 
adopted to perform our analysis is described in Section 3. Results 
are given in Sect. 4 and a summary in Sections 5. Throughout 
this work we assume a standard cosmological model with - 
0.3, Q.A = 0.7 and - 70kms"' Mpc"'. Magnitudes are given 
in the AB system. 



2. The first epoch VVDS sample 

The primary observational goal of the V VDS as well as t he sur - 
vey strategy are presented in detail in iLe Fevre et al.l ( l2005h . 
Here we stress that in order to minimise possible selection bi- 
ases, the VVDS has been conceived as a purely flux-limited sur- 
vey, i.e. no target pre-selection according to colours or compact- 
ness is implemented. In this paper we consider the deep spec- 
troscopic sample of the first epoch data in the VVDS-0226-04 
field (from now on simply VVDS-F02) that targets objects in 
the magnitude range 17.5 < Iab ^ 24. 

First-epoch spectroscopic observations in the VVDS-F02 
field were carried out with the VIMOS multi-object spectrograph 
dLe Fevre et al.ll2005l) . The LRred grism was used with 1" wide 
slits to cover the spectral range 5500 A< A < 9400 Awith an 
effective spectral resolution R ~ 230 at /I = 7500 A. The spec- 
troscopic targets were selected from the photometric catalogues 
using flie VLT- VIM OS Mask Preparation Softwai-e (VVMPS; 
iBottini et al.l l2005h . The spectroscopic multi-object exposures 
were reduced using the VIPGI tool ( Scodeggio et al. 20051 
IZanicheHi et al.ll2005l). Further details on o bservations and data 
reduction are given in lLe Fevre et al.l ( l2004aab) . 

The first-epoch data sample in the VVDS-F02 field extends 
over a sky area of 0.7 x 0.7 deg^ and has a median redshift of 
z ~ 0.76. It contains 7267 objects with secure redshifts, cor- 
responding to an average random sampling rate of 23% of the 
initial photometric sample. It is important to stress that we ex- 
plored the VVDS-F02 redshift success rate for di fferent galaxy 
types ( from star-forming to quiescent galaxies) in lFranzetti et alj 
( l2007 h obtaining a similar fraction of galaxy types in the photo- 
metric sample and in the spectroscopic sample. Thus we have 
no significant bias against special types of galaxies in our survey 
as our ability to obtain redshifts for star-forming and quiescent 
galaxies are similar 

In this work we used the [OlI]/l3727 line and the 4000 A 
break as indicators of star formation activity and stellar ages re- 
spectively. As a consequence, we limited our sample to galax- 
ies within the 0.5 < z < 1.3 redshift interval, giving us a final 
number of 4277 galaxies, after removing spectroscopically con- 
firmed broad line AGNs and stars. 



This sample provides an excellent laboratory to analyse 
spectral properties of galaxies at several redshift ranges, es- 
pecially as it is accompanied by a wealth of photometric an- 
cillary data, collecte d at several tele scopes and described in 
the fo l lowing papers: McCracken et al.i (B VRI ba nds at CFHT, 



l2003[) : lRadovichetaI 



(U-band at ESO-2.2m/WFI, 2004 ); ubvrz 



bands by the CFHT Legacy Survey project ( McCracken et al 
2007), and J and K , at ESO-NTT/SOFI (.lovino et al.. .20051 
iTemporin et al.ll2006h . 



I 



3. Methodology 

In this section we will describe in some details the methodology 
adopted to obtain measurements of [OlI]/l3727 and Dn4000 from 
VVDS spectra and to derive estimates for specific star formation 
rates (SSFRs) and stellar masses. 



3.1. Spectral measurements: [Oil] and Dn4000 

Spectral features measurements have been obtained for all galax- 
ies in our sample using the platefit-vvds software package. This 
software implements the spectral feature measurement tech- 
niques described by iTremonti et al.l 12004), but takes into ac- 
count the lower spectral resolution of our data. Here we give a 
brief outline of the fitting procedure and we refer to Lamareille 
et al. (2007) for further details. As a first step, the best fitting stel- 
lar continuum derived from a grid of stellar population synthesis 
models (Bruzual & Chariot 2003, hereafter BC03) is subtracted 
from the observed spectrum. Any remaining residual is removed 
by fitting a low-order polynomial to the continuum-subtracted 
spectrum, and then all emission lines are fitted simultaneously 
with a Gaussian profile. Finally, absorption features and spectral 
breaks are measured after having subtracted emission-lines from 
the original spectrum. 

The two spectral measurements we use in this paper are the 
equivalent width of the [OlI]/l3727 doublet (EW[OlI] from now 
onwards) and the amplitude of the 4000 A break. In particu- 
lar, we adopt the so-called narrow definition for the 4000 A 
break (hereafter Dn4000), inti-oduced by iBalogh et al.i d 19991) 
and based on the ratio of the average spectral flux density 
in the bands 4050 - 4250 A and 3750 - 3950 A ar ound flie 
break. Compared to the original definition introduced by Bruzual] 
(1983) the newly proposed index has the advantage of reducing 
the effect of dust reddening, because of the narrower portion of 
galaxy s pectrum involved in the measurement (see for details 
Kauffma nnetaLll20 03 a) . 

Comparing repeated observations available for 140 objects 
in our sample, we estimated that the typical relative uncer- 
tainty affecting our measurements is ~ 27% and ~ 10% for 
the EW[OlI] and Dn4000, respectively. Given the typical signal- 
to-noise ratio of the spectra, we estimate an average detection 
threshold for the [Oil] doublet of 8 A. Since the D„4000 ac- 
curacy can depend on galaxy redshift because some fringing 
residuals contribute to increase the noise in the VVDS spectra, 
we compared the repeated Dn4000 measurements separately for 
galaxies with z < 0.9 (where the measurement is not affected by 
fringing) and for galaxies with z > 0.9 (where our measurements 
begin to be affected by fringing). The scatter of the Dn4000 val- 
ues in the two redshift bins is always within the quoted values. 
Therefore we can safely conclude that fringing does not signifi- 
cantly affect our results. 
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3.2. Spectral classification 

To study the role played by the stellar mass in regulating the ac- 
tive phase of star formation, we adopt a parametric classification 
of galaxies based on the Dn4000 index used as an estimator of 
stellar ages ( lHamiltonlll 985; Bal ogh et alJll999l) . 

We choose a value of D,i4000 - 1 .5 to distinguish between 
spectroscopic late-type (Dn4000 < 1 .5) and spectroscopic early- 
type (Dn4000 >1.5) galaxies. 

Our choice is motivated both by galaxy evolution models and 
observations. Single stellar population (SSP) models have shown 
that in a galaxy with solar metallicity experiencing an instanta- 
neous burst of star formation, the spectrum of the underlying 
stellar population reaches values of Dn4000 larger than 1.5 in 
~ 1 Gyr (Kauffmann et al. 2003a). Different results obtained us- 
ing differ ent stellar libraries are within the error of the D n4000 
measure dKauffmann et alj|2003at iLe Borgne et"ani2006l) . while 
there is a possible dependency on metalhcity for old stellar ages 
(> 1 Gyr). Galaxies with values of metallicity as extreme as 0.05 
(2.5) Z© reac h values of Dn4000 larger than 1.5 in 2.5 (0 8) Gyr 
after a burst dKauffmann et al.ll2003al) . iLe Borgne et alj (|2006|) 
investigate the Dn4000 evolution using other modes of star for- 
mation history. They show that the Dn4000 index reaches values 
larger than 1 .5 after 1 Gyr from the initial activity also in galax- 
ies with a period of constant star formation truncated 0.5 Gyr 
after their formation, while galaxies with recurring burst modes 
always keep the Dn4000 values below 1.5 in a period as long as 
9 Gyr. 

Observationally in the present-day galaxies, values of Dn4000 
smaller (larger) than 1 .5 have identified young (old) stellar pop- 
ulations with a separati on between the two populations at 1 Gyr 
dKauffmann et al.ll2003 b). This limit has been widely used in 
other observational studies (Miller & Owen 2002; Mignoli et al] 
l2005h . 

We conclude that the Dn4000, despite the uncertainties dis- 
cussed above, is still completely adequate to provide a broad sep- 
aration of the galaxy population into active and quiescent galax- 
ies, although it might not provide a very precise age estimate for 
the stellar population. 

Hereafter we refer to our spectroscopic early- and spectro- 
scopic late-type galaxies simply as early- and late-type galaxies, 
although the classification of a galaxy on colour, morphological, 
spectroscopic based schemes should be kept distinct. The dif- 
ferent terminologies should not be mixed up, despite the robust 
trends existing among morphologies, colours, gas content, star 
formation ac tivity, and other properties (for a recent review see 
lRenzinil l2006). The category of early-type galaxies, by virtue of 
their relatively simple star formation history, are generally used 
as the preferential category of objects to test galaxy formation 
and evolution and therefore one should be careful in the defini- 
tion chosen to select them (see .Franzetti et al.,2007.) . 

3.3. Star Formation Rate Estimates 

Adopting the formula proposed by 'G uzman et al.l d 19971) we can 
translate the EW[OlI] in star formation rate estimates: 

SFR(Moyr"') « 10""-^"*''^'"«"'^«°*EW[OII] , (1) 

where we used the absolute magnitudes Mb as estimated in 
lllbert et akl d2005h . Although the average accuracy of the abso- 
lute spectrophotometric calibration of the VVDS sample is better 
than ~ 10% rm.s. in optimal weather conditions, we favour the 
use of the EW[OlI] instead of line fluxes to estimate the SFR. 
In this way we minimise the dependency of our estimates on 



less optimal weathe r conditions and slit losses (see for details 
ILe Fevre et ani2005h . 

There are three main caveats for the reliability of the SFR 
conversion obtained from EW[OlI] using Eq. 1 over the whole 
redshift interval (0.5 < z < 1.3): the uncertainty on the fraction 
of active galactic nuclei (AGN) present in our sample, the cosmic 
evolution of metallicity, and the amount of extinction. 

One risk in interpreting line measurements as proxy for 
star-formation activity is the contamination arising from the 
presence of an AGN in the central region of a galaxy. In this 
case, the emission lines in the spectrum would be a mixture of 
contributions from both the AGN and star formation. We re- 
moved from our final sample the spectroscopically confirmed 
broad-line AGNs. We estimate a contamination from type II 
AGN-dominated galaxies on the order of 7% over the red- 
shift range explored. This fraction has been estimated using 
the standard line ratio diagnostic tools, i.e. [OlII]/l5007/Hj6 vs. 
[OlI]/i3727/Hj6 in the redshift interval 0.5 < z < 0.8 and new di- 
agnostics using [OlI]/l3727, [NeIII]/l3869 and H6 in the redshift 
interval 0.9 < z < 1.3 (for details we refer to Lamareille et al. 
in prep, and Perez-Montero et al. in prep.). However, given the 
very significant uncertainty that is affecting the identification of 
these narrow-line AGNs and their relatively small number, we 
have decided not to remove them from the sample. 

Another possible limitation we face in converting EW[OlI] 
into star formation estimates is related to the metallicity 
and its cosmic evolution, that could result into a systematic 
change of line intensities with redshift even if the star forma- 
tion activity were to remain constant. iMouhcine et al.l d2005h 
have shown that a variation of the metallicity by a factor 
of 4 can produce a variation of a factor of about 3 in the 
reddening-corrected [ Oil] A3727 /Ha flux ratio (see Fig. 12 in 
iMouhcine et alj|2005h . This relation however has an opposite 
behaviour for different metal regimes. We have the scenario 
in which the [OII]/13727/Hq' flux ratio increases with metal- 
licity in metal-poor galaxies which are typically low-mass ob- 
jects (log(M/Mo) < 10, ~ 0.2Zo), and it decreases with metallic- 
ity in metal-rich galaxies which are typically massive galaxies 
(log(M/Mo) > 11, 1.4Zo), see e.g. Gallazzi et al. (2005). As for 
the direct effect of metallicity evolution on line strength, there 
is still considerable debate a bout the details of the evolution of 
the mass-metalHcity relation. lSavaglio et al.1 d2005h find very lit- 
tle change in metallicity for massive, already high-metallicity 
systems, up to z ~ 1 . They find a significant overall increase of 
metallicity f or low-mass, gen erally low-metallicity objects. On 
the contrary. i Lamareille et al.l d2006) observe a general decrease 
of metallicity which is essentially independent from galaxy 
mass, while the VVDS data (Lamareille et al. in prep.) show 
some indication of a smaller change of metallicity for low-mass 
systems than for high-mass ones. Overall these studies are con- 
cordant in finding an evolution of the metallicity within a fac- 
tor of 2 up to z ~ 1 . This would imply an uncertainty of less 
than a factor of two in our star formation rate estimates. As the 
observed variations in star formation activity are significantly 
larger than this uncertainty, and the metallicity estimates we can 
obtain are also highly uncertain because of the low spectral res- 
olution and low S/N of our spectra, we prefer not to correct our 
[On] measurements for metallicity variations. 

Finally, we also have to consider the effects of dust extinc- 
tion on our SFR estimates, although they are reduced by the 
use of equivalent widths instead of line fluxes. The observed 
line and continuum fluxes involved in the determination of the 
equivalent width are both affected by the diffuse dust extinction 
within the galaxy, and the only differential effect would be due 
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to localised extinction within HII regions, that should account 
for only a minor fraction of the total dust extinction within a 
galaxy (Kennicutt 1989; Chariot & Longhetti 2001). To account 
for the dust attenuation we corrected the EW[OlI] by the amount 
of the extinction provided by our t emplate fitting using PEGASE 
jFioc & Rocca-Volmerangd 19971 see next section). The correc- 
tion values obtained with this method vary within the interval 
EW[OlI]co,T /EW[OlI] e [1 - 2.7]. 

To the caveats we already mentioned, for sake of complete- 
ness we should also add the possible uncertainties derived from 
the unknown physical condition of the ionizing gas and those re- 
lated to the difference between our sample and the one o riginally 
used to calibrate Hq;/[OII] in Eq. 1 (lGuzmanetal.ll 19971) . 

3.4. Stellar Mass Estimates 

We obtain the stellar masses for our galaxy sample by fitting 
the photometric and spectroscopic data with a grid of stel- 
lar population synthesis models generated by the PEGASE2 
population synthe sis code (Fioc & Rocca-Volmerange 1999, 
astro-ph/9912179; iFioc & Rocca-Vofm erange 1997), and using 
the GOSSIP Spectral Energy Distribution tool (Franzetti 2005). 
The mod els were prod uced assuming a Salpeter initial mass 
function ( ISalpeterlll955[) to allow an easy compa r ison with pre- 
vious analyses (e.g. ^ Brinchmann & Elhsl l2000t : iFontana et al.l 
l2004i) . The use of other prescriptions for initial mass functions 
(IMFs) does not introduce any significant difference in the red- 
shift evolution of the mass estimate in the z interval explored as 
the ratio between masses derived with e.g. a Salpeter IMF and 
the ones derived with a Chabrier IMF is close to a constant for 
most star formation histories (see Pozzetti et al. 2007, hereafter 
P07). We use a set of delayed expone ntial star formation his- 
tories (see Eq. 3 in iGavazzi et al.ll200 2') with galaxy ages, f, in 
the range from 0.1 to 15 Gyr, and star formation time-scales, t, 
between 0. 1 and 25 Gyr. Internal dust extinction is handled self- 
consistently with the star formation activity by the PEGASE2 
code and no burst component has been added. 

Stellar masses obtained with this method were compared 
with other estimates based on the fitting of the photomet- 
ric data only and pho tometric plus spectroscopic data with 
iBruzual & Charloa (l2003i) models derived using different fam- 
ilies of star formation histories. We find the various estimates 
to be consistent with each other in the redshift interval inves- 
tigated. We estimate a statistical uncertainty on the mass esti- 
mate of approximately 0.2 dex. In particular, our methodology 
is equivalent to the smooth SFHs method described in P07 even 
if t hey adopt an expone ntially decreasing star formation history 
and Bru zual & Chariot! ("2003) models, rather than a delayed ex- 
ponential star formation history and PEGASE2 population syn- 
thesis models as in this work. The photometry coverage was im- 
proved with respect to P 07, and it now i ncludes UKIDSS NIR 
(JK bands) photometry dLawrence et all 12007) and SPITZER- 
IRA C (3.6,4.5,5.8, and 8 ix m) photometry from the SWIRE sur- 
vey dLonsdale et al.ll2003h . Even if our stellar masses are well 
consistent with the ones by P07, the new photometry improves 
the mass determination avoiding the need for any further sta- 
tistical correction, as made in P07 because of the lack of NIR 
photometry in about half of the sample. 

Stellar mass completeness limits have been estimated using 
the same stellar population synthesis models described above to 
predict apparent I-band magnitudes for galaxies with different 
star formation histories (and therefore different mass to light ra- 
tios) as a function of redshift. Comparing these predicted magni- 
tudes with the magnitude limit used to select the survey spectro- 



scopic sample we estimate that we can observe galaxies above 
log(M*) ^ 9.6,9.9, 10.5 Mo forz e [0.5-0.7], [0.7-1], [1-1.3], 
respectively. These stellar mass limits correspond to a complete- 
ness > 80% for the early-type population and more than 95% for 
late-type galaxies. 

We refer to P07 for the method adopted to construct the 
stellar mass function. Here we briefly describe the most im- 
portant points. We derive the mass functions using th e clas sical 
non-pa rametric l/Vmay formalism following Schmidt d 19681) and 
iFeltenI ([1976). We use for the fit the Schechter function limiting 
ourselves above the 80% completeness stellar mass limits for the 
two different population examined. In the estimate of the mass 
function we take into account also the incompleteness resulting 
from non-targeted sources in our spectroscopic observations and 
from spectroscopic failures. 

4. Results 

4.1. Evolution of Stellar Ages with Stellar Mass 

In the local Universe there is a well established observational 
evidence that while young galaxies are preferentially low-mass 
sys tems, old galaxie s are m ostly massive systems. In particu- 
lar, iKauffmann et alj (l2003b ) identify in their Fig. 2 the popu- 
lation of the first peak of the Dn4000 distribution observed at 
Dn4000 ~ 1.3 as emission line objects with a mean stellar age of 
1-3 Gyr and low stellar mass (log(M/Mo) < 10.4). The second 
peak at Dn4000 ~ 1.85 corresponds instead to old (> 10 Gyr), 
massive ellipticals. 

Figure 1 shows the histogram of the Dn4000 distribution in 
different intervals of stellar masses and redshifts as indicated on 
the top and on the right, respectively. For the VVDS-F02 data 
(rows 2-4) the histograms are corrected for non-targeted sources 
in our spectroscopic observations and for spectroscopic failures 
(the corrected total number and the actual number of observed 
galaxies are reported in the first and second lines, respectively, 
at the top of each panel). 

To allow a direct comparison with the local Universe, the 
first row of Fig. 1 shows the Dn4000 di stribution as derived from 
a local sample of SDSS DR4 galaxies ( Adelman-McCarthv et al.l 
2006) within the same interval of stellar masses as ours (properly 
rescaled to account for the different IMFs). This plot enables us 
to explore the presence of any relation between stellar ages and 
stellar masses in the distant Universe. The vertical dotted lines 
represent the cutoff at Dn4000 = 1.5 adopted to sub-divide our 
objects in early- and late-types, as described in the Sect. 13. 21 The 
vertical dashed line in each panel is the median of the Dn4000 
distribution in that particular stellar mass interval and redshift 
bin (see Tab. 1). In the low-mass interval at high redshift the in- 
completeness for early-type galaxies affects the high value tail 
of the Dn4000 distribution, while it does not strongly affect the 
median values, dominated by late-type galaxies. 

The first evidence shown by Fig.[T]is the different behaviour 
of the Dn4000 distribution for the five bins of stellar mass ex- 
plored. They show a clear dependency of the stellar age as in- 
ferred by Dn4000 values on the stellar mass that holds at all 
epochs up to z ~ 1.3. Even at the highest redshift bin inves- 
tigated the Dn4000 distributions closely follow the behaviour 
observed in the local Universe. Low-mass galaxies have small 
Dn4000 values and, as mass increases, the galaxy distribution 
moves to larger Dn4000 values. At any fixed redshift bin, the rel- 
ative abundance of early-type galaxies (Dn4000 > 1.5) increases 
as mass increases. The largest Dn4000 values, and thus the oldest 
underlying stellar populations or early-type galaxies following 
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Fig. 1. Histograms showing the fraction of VVDS-F02 galaxies as a function of Dn4000 in five different ranges of stellar mass and in 
three redshift bins (0.5-0.7, 0.7 - 1 .0, 1 .0- 1 .3, r ows 2-4). We plot as a compariso n in the first row the Dn4000 galaxy distribution for 
the local Universe (z ~ 0) from the SDSS DR4 jAdelman-McCarthv et alj|2006l) . The counts correct ed for non-targeted sources in 
our spectroscopic observations and for spectroscopic failures (unidentified sources) following lllbert e t al. (2005) are reported in the 
first line at the top of each panel, the observed counts of (total, late-, and early-type) galaxies in the second line. The vertical dotted 
lines at Dn4000 =1.5 represent the adopted separation between early- and late-type galaxies. The vertical dashed lines represent 
the median Dn4000 distribution in each panel. The panels in the first column dominated mostly entirely by late-type galaxies have 
dashed Dn4000 distribution when the mass completeness for the late-type population is less than 50%. 



Sect. |3.2| are hosted in the most massive galaxies at aU redshifts 
up to z ~ 1.3. 

Another striking evidence coming from Fig.[T]is that for the 
lower mass galaxies (log(M/Mo) < 10) the shape of the Dn4000 



distribution does not seem to evolve with time between z ~ I 
and z ~ 0- For these galaxies the Dn4000 distribution is peaked 
around Dn4000 ~ 1 . 1 in all redshift bins, suggesting a good cor- 
respondence up to high redshifts between a low stellar mass and 



6 



D. Vergani et al.: Galaxy assembly 



Table 1. Dn4000 median distribution 



lO'-'-^ Mo lO'-S-'" Mo lO'"-'"^ Mo 10'"^-" Mo 10^" Mo 

z 



~0 1.24 ±0.00 1.28 ±0.00 1.46 ±0.00 1.66 ± 0.00 1.81 ± 0.00 

0.5 - 0.7 1.18 ± 0.01 1.20 ± 0.01 1.27 ± 0.01 1.50 ± 0.02 1.69 ± 0.02 

0.7- 1.0 1.13 ±0.01 1.16 ±0.01 1.21 ±0.01 1.31 ±0.01 1.55 ±0.02 

1.0 - 1.3 1.07 ± 0.04 1.07 ± 0.02 1.09 ± 0.01 1.20 ± 0.02 1.37 ± 0.02 



NOTES: Dn4000 median values in five intervals of stellar mass 
and three redshift bins for our VVDS-F02 sample. The error 
quoted in each interval represents the error on the median dis- 
tribution using a statistical Jackknife technique. The median val- 
ues of the Dn4000 distribution in a nearby galaxy sample (SDSS 
DR4) are also Usted. 



a young underlying stellar population. Intermediate and high- 
mass galaxies show significant evolution with redshift at z > 0.7 
as testified by the progressive shift of the median value (dashed 
line in Fig. 1). As cosmic time goes by, for stellar mass larger 
than log(M/MQ) > 10 galaxies populate progressively the locus 
of larger Dn4000 values and a secondary peak emerges in the 
Dn4000 distribution. The gradual accumulation of the secondary 
peak is a process started at early epochs that continues efficiently 
down to the local Universe. 

We find an age-stellar mass relation with low-mass galaxies 
having younger stellar populations and more massive galaxies 
possessing older stellar population. This correlation, originall y 
observed in the local Universe by e.g. lKauffmann et alj (l2003bl) . 
seems to hold up to z - 1.3. Our findings on the progressive 
accumulation with time of massive galaxies with large Dn4000 
values can be interpreted in the framework of a top-down picture 
of the stellar mass assembly history of galaxies, also termed as- 
sembly downsizing. With this terminology we refer to the galac- 
tic property that regardless of type, massive galaxies are fully 
assembled at earlier times, with less massive galaxies assembled 
later (when mass is defined as stellar mass). Note that here we 
make no statement about how the stellar mass was aggregated, as 
both merging and star formation can contribute to the assembly 
process. However we will argue in Sect. l4.5l that galaxy major 
merger does not seem to be a fundamental mechanism, at least 
atz < 1. 

Evidence of this top-d own stellar ma ss assem bly is already 
presen t in hterature , e.g. Thomas et al. (2005); Cima tti et alj 
( l2006l) : iBundv et all ( 12006 ). Studies on stellar age, metallicity, 
and cosmic star formation at z ~ prove that stars in massive 
galax ies were formed long ago and over a short period of time 
(e.g. iThomas et alj 12001 iHeavens et all 12004 iJimenez et al] 
l2005h . At high redshifts the reddest galaxies already possess 
old a ges (> 1 Gyr), establishing their forma tion epoch at 
z > 2 dMcCarthv et alj|2004l: Icimatti et alj|2004l) . Furthermore, 
massive proto-disk galaxies already ex ist at z ~ 2-3 (e.g. 
Labbe et al. 2003; Stockton et al. 2004; F orster Schreiber et alj 
20061; iGenzel et alj|200 6). On the contrary, the underlying stel- 
lar ages of late-type galaxies apparently do not show evolution. 
This observation is probably justified by a continuous and effi- 
cient formation rate of new stars (see next Sect. l4.3b as discussed 
also in other studies ( Arnouts et al. 2007; Faber et al. 2007). 

Theoretical models, however, predict that stars in an object 
can be born in a different place from where they finally end up 



(e.g. de Lucia et al. 2006). Even if we observe a downsizing sig- 
nal, stars born long ago in several objects could still coalesce into 
a single object in a hierarchical way making the concepts of as- 
sembly of a system and assembly of its stars different. However, 
to be able to test this theoretical prediction we need to wait for 
the next generation of surveys and more detailed semi-analytical 
models. 

Based on our phenomenological approach we can conclude 
as follows. Our sample shows to which extent the stellar mass 
drives the galaxy evolution. For the first time we present a direct 
comparison of stellar ages and stellar masses of galaxies over 
the last 8 Gyr ex tending to higher redshift the seminal work by 
.Kauffmann et a n f2003b.) on z = galaxies. 

4.2. SFR dependence on stellar mass and redshift 

A cosmic decline of the star formation activity has been pre- 
sented in several st udies, startin g from field redshift surveys by 
iLilly et al.1 (Il996l) : iMadau et aP (Il996h . Nowadays, the knowl- 
edge of the cosmic star formation history out to redshift z ~ 1 
is quite profound (Hopkins & Beacom 2006, and references 
therein). However many aspects of galaxy evolution, like the 
contribution to galaxy mass assembly of various categories 
of galaxies, still n eed to be understood (see for details e.g. 
ICimatti et al.l l2006). 

Panel a of Fig.|2]shows the median distribution of EW[OlI] in 
the redshift interval z = 0.5-1.3 as a function of the stellar mass. 
Galaxies with no EW[OlI] detection are taken into account in the 
median computation assigning to them the EW[OlI] upper limit 
estimated from the 1 cr measurement uncertainty. Interpreting 
the EW[OlI] as a proxy for star formation activity, we can in- 
fer that the latter has a clear redshift evolution which is strongly 
modulated by stellar mass. We observe the already well estab- 
lished global star formation decline from z ~ 1 to the present 
epoch, and also that within a given redshift bin the smaller the 
stellar mass of the galaxy, the more actively it is forming stars. 

These results are in excellent agreement with the findings by 
iFeulner et al.1 (12005 a) based on a large sample of 9000 galaxies 
although with only photometric redshifts from the FORS Deep 
Field and the GOODS-S field. They show a similar trend for 
the SSFRs at different stellar mass intervals with no contribu- 
ti on to the stellar mass growth for massive galaxies (cf. Fig. 2 
in'Feulner et al."2005a'). We also find an excellent agreement in 
the common redshift rang e with the study of 2 07 galaxies from 
the spectroscopic GDDS d Juneau et al.ll2005h . They show that 
massive GDDS galaxies (log(M/Mo) > 10.2) are in a quiescent 
mode, and lower mass GGD S galaxies in a bur st phase. Similar 
results have been obtained by Zheng et al. ( 2007) using ~ 15,000 
COMBO- 17 galaxies with photometric redshifts, supplemented 
by ultraviolet and infrared data to account for otherwise unde- 
tected star formation activity. 

Our results, being obtained on a large spectroscopic sample, 
put on firmer grounds the results previously found using photo- 
metric galaxy catalogues or smaller spectroscopic surveys: the 
star formation in massive systems is at very low values since 
z ~ 1-1.3, while star formation increases when considering less 
massive systems at each redshift bin explored. 

Panels b and c of Fig.|2]show what happens when one consid- 
ers early- and late-type galaxies separately. Very massive galax- 
ies (log(M/M0) > 1 1) irrespective of their spectral type show 
negligible star formation activity at each redshift bin probed. 
Moving to lower masses, galaxies of different spectral classifica- 
tion show different behaviour (both in slope and normalization 
of EW[OlI] values). This result holds over the entire inspected 
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Fig. 2. The median EW[OlI] as a function of the redshift for our 
VVDS-F02 sample (a) in bins of stellar masses, coded as fol- 
lows: 9.0 < log(M/M0) < 9.5 (blue, circle), 9.5 < log(M/Mo) < 
10 (cyan, square), 10 < log(M/Mo) < 10.5 (magenta, triangle), 
10.5 < log(M/M0) < 11 (red, diamond), and log(M/M0) > 11 
(green, pentagon). The same visualization is presented for the 
spectroscopic early- (b) and late-type galaxies (c). The symbols 
represent the observed EW[OlI] and the arrows the dust cor- 
rected EW[OlI]. The points are connected with dashed (dotted) 
lines when the mass completeness of each population is above 
(below) 50%. 



redshift range. In panel b at masses lower than log(M/Mo) ~ 10 
there is an indication of the presence of galaxies with an old 
underlying stellar population (ETGs) and non-negligible lev- 
els of [Oil]. However, given a number of considerations, in- 
cluding the larger mass incompleteness and the small num- 
ber statistics for this class of objects (see numbers quoted in 
Fig. 1), one should be cautious in putting too much weight on 
this result. Interestingly, with our classification scheme we find 
a lower contamination of star forming galaxies among early- 
type objects than reported in previous studies (5% vrs. 20-50%, 
see Strateva et al. (200jt; iFranzetti et all (l2007l) : ICimatti et all 
( 12002) : .Kriek et al.i (12006*)). This could be the result of our clas- 
sification criterion that, being based on Dn4000, allows reject- 
ing e.g. dusty-star forming EROs that p ossess red colours but 
low D„4000 values (MignoU et aLll2005h . Another concern we 
need to stress is the possible contamination from the type II 
AGN-dominated galaxies as they might still be present in our 
sample. A number of studies have recently shown that emission 
lines in early-type galaxies at z - are mainly due to AGN 
rather than to star formation activity (e.g. IWeiner et al.ll2007l) . 
Therefore we can not exclude the possibility that even the 5% 
of early-type galaxies with non-negligible star formation could 
actually be AGNs, thus casting doubt on the possibility of inter- 
preting EW[OlI] as star formation indicator for low-mass early- 
type galaxies in panel b of Fig. 2. 



Panel c of Fig. 2 shows clear trends for late-type galaxies 
both as a function of redshift and of mass. Galaxies have a 
stronger EW[OlI] going back in cosmic time and this trend is 
steeper for less massive galaxies. The trends visible in panel a of 
Fig. 2 are therefore driven by the late-type galaxies. 

An important piece of information about the global star for- 
mation history of galaxies becomes apparent when comparing 
the results summarised in Figs.[T] and |2] It is known that the 
Dn4000 strength is the result of the cumulative star formation 
history of a galaxy. As already discussed in Sect. 13. 21 Dn4000 
values below 1.5 are indicative of significant star formation ac- 
tivity within the last 1 Gyr, while values above 1.5 are associ- 
ated to galaxies that have finished forming their stars at least 
1 Gyr ago, and are therefore evolving mostly, if not exclusively, 
via the passive evolution of their stellar population. Instead the 
[Oil] line is sampling only the instantaneous star formation ac- 
tivity, because the ionizing flux in HII regions is provided almost 
entirely by massive stars (log(M/M0) > 10) with lifetimes lim- 
ited to 20 Myr or less. Therefore, the star formation time-scales 
probed with the Dn4000 and the [Oil] equivalent width are very 
different. Despite their diff'erence, however, there is a general 
agreement between them: both indicators consistently point to- 
wards a star formation history peaked in the past for massive 
galaxies, with little or no residual star formation taking place at 
the epoch of the observation. Vice-versa most of the low-mass 
systems show just the opposite characteristics: almost all galax- 
ies in this subset have small Dn4000 and a rather strong current 
star formation activity. 

In our sample massive galaxies are the only subset of the 
galaxy population that shows a significant presence of objects 
with large Dn4000, and all these objects have also very low or 
null recent star formation activity. For low-mass galaxies both 
indicators suggest a general picture of extended star formation 
activity. These galaxies are dominated by the young stellar pop- 
ulation at all epochs (although with this kind of analysis we can- 
not say anything about their epoch of formation, of course). 

The very good agreement observed for massive objects be- 
tween the instantaneous and the integrated star formation indica- 
tors can be understood by assuming that the dominant mode of 
star formation activity for these systems after an initial burst is 
a continuous one, with relatively smooth variations in time, and 
with no significant contribution from a secondary stro ng burst. A 
similar conclusion has been presented recently by Noes ke et al.l 
(2007), based on their analysis of the AEGIS survey data. 

We notice that a somewhat similar dependence of burst ac- 
tivity on galaxy st ellar mass has been observed in the SDSS 
data by Kauffma nn et al.l ( 2003al) . although their most recent 
work (IKauffmann et al.l 12006 *) indicates a more important role 
for galaxy surface mass density than for total mass in regulating 
star formation activity. 

The agreement between the two star formation indicators 
would be violated by the small fraction of low-mass spectral 
early-type galaxies where a large Dn4000 value, that points to- 
wards a predominantly old stellar population, is coupled with a 
significant [Oil] emission, should this feature be indicative of 
ongoing star formation activity (see Fig. 2, panel b). We have al- 
ready listed a series of caveats on the risks of over-interpreting 
the significance of such population. It is nevertheless worth men- 
tioning that it is possible to reproduce such a dichotomy between 
the two indicators when the starburst is involving only a very 
small fraction of the galaxy total stellar population (approxi- 
mately less than 2% by mass, but the exact value for this limit 
depends on the amount of dimming introduced by dust extinc- 
tion localised in the star-forming regions). 
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The general picture which emerges from the Figs.[T] and|2] 
follows the classical concept of downsizing as firstly suggested 
by[Cowie et al. ( 1996). In this framework one expects a continu- 
ous growth in stellar mass of low-mass galaxies with time while 
the massive galaxies have stopped to form stars at early epochs. 
This galax y attitude is als o termed downsizing in star formation, 
or in time (iNeistein et al.l |2006). Most of the claims of downsiz- 
ing in the literature fall in this category ( e.g. iBauer et al. ||2005t 
lJuneau et al.|[2005l: iFeulner et al.ll2005alltt iTn 



resse et al.ll2006ir 



4.3. Stellar Mass Assembly in Early- and Late-Type Galaxies 

We show in Fig. [3] the stellar mass function (MF) in three red- 
shift bins both for early- and late-type galaxies according to the 
spectral classification introduced in Sect. l3.2l We also show the 
total mass function for each redshift bin, but we refer to P07 
for a detailed discussion on the VVDS-F02 mass function. The 
parameters of the Schechter fit (a, M*,^,^, (p*) for early- and late- 
type galaxy samples are listed in Tab. 2 along with their uncer- 
tainties. As the interval 0.7 < z < 1 provides the best sampling 
of the overall galaxy population we adopt the characteristic mass 
(M*,^,.) of late-type galaxies measured in this bin for the lower 
redshift one (where we do not sample well enough the high-mass 
end of the population), and the slopes of both early- and late-type 
MFs measured in this bin for the higher redshift one. The verti- 
cal dashed line in each panel represents our mass limit in that 
redshift bin (see Sect. 13.4b . To judge the evolution in the three 
redshift intervals we over-plot the Schechter fits to early- and 
late-type stellar mass functions of the lowest redshift bin at each 
panel (red and blue, dotted lines), and the total local MF range, 
using Cole et al. (2001) and Bell et al. (2003) (shaded region and 
dotted line, respectively). 

Comparing the total galaxy MF of Fig. 3 to the local MF, we 
find a strong evolution depending on galaxy stellar mass con- 
firming previous results of P07. Galaxies with stellar masses 
larger than log(M/M0) > 1 1 show an evolution on the order of 
30 - 50% up to z ~ 0.7-1.0 and a faster evolution above z ~ 1, 
while the number density of less massive galaxies decreases 
more continuously with redshift. 

Our main focus is on the contribution of the diff'erent galaxy 
populations to the total mass function and to its evolution. For 
stellar masses larger than log(M/Mo) > 1 1 .4 most of the galax- 
ies have an early-type classification over all the redshift range 
explored (0.5 < z < 1.3) and no massive late-type galaxies 
(log(M/M0) > 1 1.5) are observed at z < 0.7 within our sam- 
ple. Figure 3 also shows that at the intermediate/low-mass tail 
(log(M/M0) < 10) the number density of both early- and late- 
type galaxies increases with cosmic time: by a factor of ~ 1.6 at 
log(M) =10 Mq within the redshift range explored (and a fac- 
tor of 3 compared to local MFs) for late-type galaxies, and more 
strongly for early-type ones (by a factor of 5 within the redshift 
range explored). We can therefore conclude that while the mas- 
sive end of the mass function is mainly dominated by early-type 
galaxies up to z ~ 1, late-type galaxies mostly contribute to the 
intermediate/low-mass part (< 10.5 M©) of the mass function at 
all redshifts. Furthermore the abundance of massive early-type 
galaxies (> 1 1 M©) increases with cosmic time but with a de- 
clining rate at lower redshift. There is also an indication of a 
small decreasing of the massive tail of the late-type mass func- 
tion from z- 1.3 to z = 0.5. 

According to our results, the stellar mass function shows 
an increasing (decreasing) relative contribution of massive early 
(late) type galaxies, and a general increase, but at different rates 
and depending on the mass, of both early and late-type galax- 
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Fig. 3. Stellar mass function for the spectroscopically classified 
late-type (empty circles) and early-type (filled circles) galaxies 
from VVDS-F02 sample. The total mass function for each red- 
shift bin is plotted as dashed-point line, while in each panel we 
also report the total local MF range, using Cole et al. (2001) and 
Bell et al. (2003) (shaded region and dotted line, respectively). 
The stellar mass functions of the two spectral types of the lowest 
redshift bin are over-plotted at each panel as dotted lines. The 
vertical dashed lines represent the 80% completeness limits for 
the two spectral types. 



ies with cosmic time. These results imply a redistribution of the 
stellar mass amongst the spectral types with a declining contribu- 
tion to the overall mass function of massive late-type with time 
associated with an increasing number density of intermediate- 
massive early-type galaxies. 

The overall picture emerging from our stellar mass function 
agr ees quite well with previous results (e.g. iFontana et al.ll2004l : 
Bun dv et"an i2006). Generally the stellar mass function of late- 
type galaxies agrees very well in all studies. As far as the stel- 
lar mass function for the early-type galaxies is concerned the 
better agreement is obtained with the study by iFontana et al.l 
(2004) who use a spectral classification, instead of a rest-frame 
colour scheme (Bundv etal. 2006). A slightly worse concor- 
dance is found for the early-type galaxy densities at z > 0.7 with 
Bundy et al. which are larger than the one computed in our study. 
Furthermore, we find a stronger evolution of early-type massive 
galaxies at z > 1 compared to the one measured for red galaxies 
in Bundy et al. 

The mass redistribution of morpholog ical types as a func- 
tion of redshift was first noticed by iBrinchmann & Ellii 



Table 2. Schechter Fit Parameters 
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z range 




a 


{hiQ M0) 


{\Qi-^hl,.Mpc-^) 




a 


(hiQ Mo) 


(10-^ h^^„Mpc-^) 


logM,,. 


0.5 < z < 0.7 
0.7 <z< 1.0 
1.0 <z< 1.3 


ETG 


-0 36+" ''^ 
-0 46+"-^'' 
-0.46 


11 06+"'^ 
J ^ 09+0.12 

1 1 20+^-'^^ 


1.50:|j^ 

0.80:1 

0.25i!^ 


LTG 


_1 44+0.11 

_i 14+8:?^ 
^•■^ -0.12 
-1.34 


10.96 

10.96!«:'^ 


1 05+'" 
1 30+'''' 
75+'' ' 


10.87 
11.11 
11.42 







NOTES: The listed values are computed on the early-type (ETG) and late-type (LTG) samples using a classification scheme 
based on the age-estimator Dn4000 (see Sect. |3.2t . 



(l2000l) computing the stellar mass density in a small sam- 
ple of I-selected galaxies with spectroscopic redshifts and in- 
frared photometry^JJsing Hubble Space Telescope morpholo- 
gies, iBrinchmann & Ellii (2000) observed a decline in the 
mass density of iiTegular galaxies from z ^ 1 to today ei- 
ther associated to a transformation into regular morpholo- 
gies and/ or to merging mech anisms. These results were con- 
firmed by iBundv et alJ (120051) using a sample of approximately 
two thousand morpholo gically-cl assified galaxies from GOODS 
dGiavahsco et al.ll2004l) . Later on iBundv et"aLl (|2006|) provided 
further support to the stellar mass redistribution using the 
DEEP2 spectroscopic sample producing type-dependent stellar 
mass functions based on rest-frame colours. Further evidence 
can also be found in other studies. For example, a similar mod- 
est change in abundance of massive early-type up to z ~ 1 was 
found in the K20 sample which has a high redshift accuracy and 
a spectral classification similar to ours (Fontana et al. 2004). The 
constancy in blue populations and the increasing trend of the 
red sequence classified according to (NUV - r) colours have 
been observed in a study of stellar mas s density by Arnou ts et 
al. (2007), and pre viously observed by iBorch et all (l2006l) and 
[Martinet al.l(l2007l) . 

We observe in Fig.|3]a downward transition with time of the 
threshold in the stellar mass above which the early-type galaxies 
become dominant in the relative contribution to the total stel- 
lar mass function. This characteristic mass, or transition mass 
(Mfr reported in T ab. 2), was detected in the local Univers e 
around 3x10 ' " Mm (iKauffmann et al.l2003bl;lBaldrv et al.l2004l) . 
iBundv et al.l (l2006 l) report its evolution with cosmic time us- 
ing three different methods to define the two galaxy categories: 
galaxies are partitioned according to rest-frame colours, to the 
star formation rate and to morphological properties. Our findings 
based on the Dn4000 index are in between those ob tained divid- 
ing the galaxies by colours and by morphology in i Bundv et al.l 
(I2OO6I) and show a similar trend in cosmic time. A less steep evo- 
lut ion is observed wh en considering the classification obtained 
bv lBundv etsl\ ( |2006|) with the SFR indicator. 

To better show the differential evolution of galaxies and the 
general redistribution in stellar mass among galaxies we plot 
in Fig.|4] the evolution of the early-type galaxy fraction over 
the entire (early- and late-type) population for different stellar 
mass bins as a function of redshift. Both observed counts and 
Schechter fitted values are used in the computation and plotted 
with different symbols. The eiTor-bars include both the uncer- 
tainties in the spectral measurements and the effect of changing 
the classification due to the D,i4000 errors. As in Fig.[T]we report 
the reference points at z ~ from a local galaxy sample (SDSS 
DR4). 

Figure 4 clearly shows how the evolution of the relative 
abundances of early- and late-type galaxies strongly depends 



on galaxy stellar mass. The relative fraction of the most mas- 
sive galaxies (log(M) > 11 M©) significantly changes with cos- 
mic time, between z = 1.2 and z = 0.85 due to both the de- 
crease of massive late-type and the increase of massive early- 
type galaxies, followed by a flattening below z = 0.85, when 
this second population becomes dominant in this mass regime. 
A similar flattening is present at z < 0.6 for galaxies with 
log(M) > 10.5-11 Mq. a more continuous evolutionary trend 
is suggested for the relative fraction within the less massive pop- 
ulation (< 10.5 Mq): both late- and early-type densities increase 
but the first category of galaxies is always dominant at this stellar 
mass regime (Fig.|3ll. 
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Fig. 4. The evolution of the fraction of galaxies classified spec- 
troscopically as early-type in different mass intervals. The mass 
bins are coded as follows: 9.0 < log(M/M0) < 9.5 (blue, circle), 
9.5 < log(M/Ma) < 10 (cyan, square), 10 < log(M/M0) < 10.5 
(magenta, triangle), 10.5 < log(M/M0) < 11 (red, diamond), 
log(M/M0) > 11 (green, pentagon). Filled symbols represent 
the counts of observed galaxies, while the open symbols repre- 
sent the counts derived from lA'mo.v formalism. The dashed lines 
are the Schechter fitted values to the stellar mass function (dotted 
in the case of mass completeness for early-type galaxies below 
50%.) The points at z ~ are our SDSS reference points. 
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As these measured ratios could be affected both by differ- 
ent mass completeness limits for early and late-type galaxies (as 
the VVDS is a magnitude limited survey) and by any incom- 
pleteness in the global spectroscopic sample, we derive the rel- 
ative fraction of early-type galaxies also integrating the best fit 
Schechter function. The agreement between the observed ratios 
and the ones derived by the Schechter integration shows that our 
results are robust. 

To summarise, our results suggest a galaxy mass assembly 
history with a mild evolution with redshift in number density 
(P07) for very massive galaxies (> 11.4 Mq), which are mainly 
constituted by early-type galaxies up to the redshift sampled 
(z < 1.3). Since the relative fraction of massive early- and late- 
type galaxies evolves with cosmic time we witness a spectral 
transformation of late-type systems into old massive galaxies at 
lower redshift. In Sect. l4.5l we analyse the efficiency of this mass 
assembly process driven by the star formation rate observed in 
galaxies. 

4.4. Influence of cosmic variance 

An intrinsic source of uncertainty in any estimate of number den- 
sities from a finite-size survey is represented by the variance in- 
duced by density fluctuations on scales comparable to, or larger 
than, that of the sampled volume. This is what is commonly 
termed "cosmic variance". It inversely depends on the sample 
volume and directly on the two-point correlation function of the 
class of galaxies under analysis. This cosmic variance will in- 
troduce a systematic effect on the overall amplitude of our mass 
function. The significance of this effect needs to be evaluated. 
This can be obtained in our case with fairly good accuracy for 
the three redshift bins of VVDS-F02 used here, since we akeady 
know the clustering properties of both the early- and late-type 
galaxy population ov er the whole range covered by the survey 

jMeneux et aLll2006h. 

Following IPeeblesI (Il980l) and ISomerville et al.1 (|2004|) . un- 
der the assumption of a power-law correlation function ^(r) - 
(f/foy^ the relative variance of galaxy counts over a given vol- 
ume with radius R can be expressed as 

cri^JiiR/ror ■ (2) 

Here J2 is a classic al integr a l over the two-point correlation func- 
tion as defined in iPeeblesI (I1980D . which in the power-law ap- 
proximation can be expressed as 

72 =72/[(3-r)(4-r)(6-r)2n . (3) 

In the case of a non-spherical volume V, we can nevertheless de- 
fine an equivalent radius as R = (3y/4:7r)'^^. For our three sam- 
ples, this corresp onds to R - 34 , 45 and 49 h ' Mpc, respec- 
tively. From Meneux et al.l (12006 ). we know that, to sufficient 
approximation, the correlation length ro and slope 7 of ^(r) do 
not evolve significantly, being ro - 4 (2.5) h"' Mpc and y = 1.8 
(1.5) for early- (late-) type galaxies respectively, over the full 
redshift range under investigation. Thus, we have all the ingre- 
dients to use eqs. |2]and|3]to estimate the expected nns relative 

error on the measured densities, crcounts - -^J^^ for each of the 
three samples. The results are reported in Table |3] The listed 
values indicate that the systematic error on our estimates of the 
mass due to cosmic variance is expected to be between 13 and 
20%. This is fully consistent with the result we have obtained 
by splitting the sample into two parts and recomputing the mass 



Table 3. Relative rms error on galaxy number densities expected 
from cosmic variance in the three redshift samples 







Type 


0.5 <z< 0.7 


0.7 < z < 1 


1 <z< 1.3 


Early 


0.20 


0.15 


0.14 


Late 


0.17 


0.14 


0.13 



functions separately for the two subsamples. The agreement be- 
tween these two independent estimates allows us to conclude 
that our conclusions are robust with respect to cosmic variance. 

4.5. Efficiency in the Stellar Mass Assembly Process 

In the case of dry merging - merging of quiescent galaxies with- 
out gas (star formation) involved - the age of the galaxy (i.e. the 
time elapsed since it assembled most of its current stellar mass) 
may not coincide with the age of its underlying stellar popula- 
tion. As a consequence the processes of star formation and mass 
assembly remain distinct: it is possible that stars which formed 
long ago were assembled only recently in a newly crea ted galaxy 
( lRenzinill2007Hcimatti et al.ll2006l: iBundv et al.ll2006l) . 

Within such a scenario one considers two separate downsiz- 
ing signatures, i.e. downsizing in star formation (the transfer of 
star formation activity to lower mass galaxies) and downsizing 
in mass assembly (irrespectively by type classification, massive 
galaxies are fully assembled at early times, with less massive 
galaxies assembling later). If one or more mechanisms are re- 
sponsible for quenching star formation activity earlier in massive 
galaxies, and if dry mergers represent only a minor contribution 
to the galaxy mass assembly, then the two concepts of down- 
sizing express one single idea. In the previous sections we have 
found that the distribution of the stellar ages and stellar mass 
can be interpreted in the framework of a top-down picture of the 
stellar mass assembly history of galaxies. However, we made no 
hypothesis on the mechanism or mechanisms that lead the stellar 
mass to aggregate, leaving open the possibility of mass accretion 
by both merging and star formation activity. 

Can the galaxy star formation rates and their efficiencies 
alone justify the observed mass assembly history without invok- 
ing any (dry) merger mechanism? To investigate this issue we 
define for each galaxy at redshift z the time At as the time in- 
terval between the age of the Universe at the observed redshift 
z and the time corresponding to the lower limit of the redshift 
bin in which the galaxy falls. For three redshift bins in Fig. 5 we 
plot with left-leaning diagonals (blue-coded) the Dn4000 distri- 
bution of galaxies which assemble enough stellar mass to move 
to a higher mass interval only on the basis of their estimated 
SFRs within a time At. The assumption is that galaxies will sus- 
tain a SFR, equal to the observed one, over the time At. This 
appears to be a reasonable assumption as our At is never larger 
than 1.3 Gyr. The right-leaning diagonals (red-coded) represent 
the galaxies which, given their observed SFR, would remain over 
the time Af in the mass interval assigned at the time of the ob- 
servations. 

Figure |5] shows that the lower the stellar mass the larger is 
the fraction of galaxies which sustain a significant stellar mass 
growth, i.e. not only forming stars as a local and episodic event 
but also effectively assembling mass. Once a galaxy reaches a 
certain stellar mass, it terminates to grow efficiently. The phys- 
ical reason could be the natural exhaustion of the gas reservoir 
contained in massive galaxies which did consume it efficiently in 
the first phase of their life. On the contrary, low-mass galaxies. 
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9.0<log(M)<9.5 9.5<log(M)< 1 1 0<log(M)< 1 0.5 1 0.5<log(M)< 1 1 log(M)>11 




0.5 1 1.5 20.5 1 1.5 20.5 1 1.5 20.5 1 1.5 20.5 1 1.5 2 

D„(4000) 



Fig. 5. Histograms showing the fraction of VVDS galaxies as a function of Dn4000 for five different ranges of stellar mass at three 
redshift epochs (0.5 - 0.7, 0.7 - 1 .0, 1 .0 - 1 .3). The corresponding look-back time is also indicated. We define Af for each galaxy 
at redshift z as the time interval between the age of the Universe at the observed redshift z and the time corresponding to the lower 
limit of the redshift bin in which the galaxy falls. Left-leaning diagonal (blue) histograms represent those galaxies that assemble 
within a time Af a stellar mass enough to move to a higher mass interval and with right-leaning diagonals (red) those staying in the 
mass interval assigned at the time of the observations. The assumption is that these galaxies will sustain a constant observed SFR 
over the time Af, which turns to be a reasonable assumption as our Af is never larger than 1.3 Gyr. 



having sustained a lower rate of star formation activity over the 
time still have a certain amount of fueling. 

The fraction of galaxies that increase significantly their stel- 
lar mass (blue-coded population) is larger than the fraction of 
galaxies that is not effectively assembling mass (red-coded) at 
lower stellar masses with time. The evolution of the character- 
istic stellar mass at which blue- and red-coded populations are 



equally contributing to the total population parallels that of the 
quenching mass, defined as the mass at which n o active systems 
are observed (see for details lBundv et alj|2006l) . 

For each stellar mass bin Fig. 5 shows a decline in the mass 
assembly efficiency from high to low redshifts. Similarly for 
each redshift bin there is a decline in the mass assembly effi- 
ciency moving from low to high stellar masses. Interestingly, the 
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few massive galaxies which can still potentially increase their 
mass are classified as late-types from Dn4000. 

These results and those obtained in previous sections sup- 
port the hypothesis of apparent no evolution in the blue se- 
quence proposed by Faber et al. (2007) and the mild evol ution of 
the ea rly-type class of galaxies (Zucca et al. 2006; Brow n et"al] 
l2007l Arnouts et al. 2007) which are largely inefficient at as- 
sembling mass at all epochs (at least since z ~ 1.3), and are pro- 
gressively decreasing their efficiency as cosmic time progresses. 
In particular, on the basis of the mass migration from blue to red 
objects, Arnouts et al. (2007) propose a similar time-scale for the 
period required for blue galaxies to quench their star formation 
activity and move to the red sequence, and the birth rate time of 
blue galaxies which are continuously refilled by new stars. 

If we consider the galaxies with stellar mass M, 
(log(M/M0),=i,2,3,4,5 = 9-9.5, 9.5 - 10, 10- 10.5, 10.5 - 1 1 ,> 1 1) 
observed at redshift zj (Zj=i,2,3 = 0.5 - 0.7, 0.7 - 1 , 1 - 1 .3), their 
progenitors are the sum of two galaxy contributions if no merg- 
ing activity is advocated. The first progenitor class is constituted 
by galaxies with mass M,_i at zj-i that are assembling enough 
stellar mass to migrate in the higher mass interval within a time 
Af, where At was previously defined (e.g. blue-coded galaxies at 
z - 0.7-1 and log(M) - 10.5-11 Mq). The second contributors 
are galaxies with stellar mass included in the same mass interval 
(M,) at the epoch zj-i that in a time At are not growing suffi- 
cient mass to flow in the subsequent stellar mass interval (e.g. 
red-coded galaxies at z = 0.5-0.7 and log(M) > 1 1 Mq). 

We compute within our mass completeness limit the number 
(per unit of co-moving volume) of galaxies observed at redshift 
Zj of stellar mass M, and that of progenitors at redshift zj-i as de- 
scribed above. We obtain that the number of progenitors at Zj-\ 
defined as above can account, on average, for (80+10)% of the 
galaxies at Zj of stellar mass M,, and therefore that no signifi- 
cant contribution to the mass assembly process is required from 
merging activity. The agreement between number of progeni- 
tors and number of observed galaxies improves moving to higher 
masses, reaching nearly 100% forlog(M) > 10.5 Mq. This result 
shows that star formation activity alone can explain the observed 
evolution of stellar mass function up t o z ~ 1 . A similar con- 
clusion has been presented recently by i Bundv et al.l (|2007|) for 
spheroidal galaxies comparing the dynamical and stellar mass in 
the GOODS fields. The match between the growth rate in the 
abundance of spheroidals with that predicted by the assembly of 
dark matter halos enables Bundy et al. to state that merging is not 
the primary mechanism to form spheroids and additional mecha- 
nisms as morphological transformations are required to drive the 
observed evolution. Further convincing evidence that a signifi- 
cant number of massive blue galaxies must have quenched their 
star formation and moved to the red seque nce without invoking 
any merging mechani sm is provided by e.g. lScarlata et al.l(l2007h 
and lBell et al.l (l2007l) . This picture satisfies both the luminosity 
and mass function constraints. 

The scenario described here, if correct, implies that all recent 
observations on galaxy evolution can be summarised with one 
single concept of downsizing, as downsizing in star formation 
and downsizing in mass assembly go in parallel. Still, a more 
detailed and quantitative analysis of the mass function and star 
formation rate evolution are needed before a more accurate eval- 
uation of the role played by mergers in the build up of the galaxy 
mass can be obtained. It is in fact well known that mergers do 
take place and recently individual ca ses o f dry mergers have also 
been Hsted by e.g. Be U et al.1 (|2006|) and lTran et all (|2005|) . We 
cannot therefore exclude the possibility that between two red- 



shift bins a small fraction of galaxies moves from mass bin M,- 
to a higher mass bin via merger. 

A possible objection to the scenario in which merging would 
not be important is the large fraction of galaxies up to z = 1 
possessing spatial and kinematic asymmetries (e.g. lFlores et al.l 
2006; Puech et al. 2007). These lopsided galaxies are very often 
associated to merging of galaxies. However there exist various 
ranges of mass merging fractions and several mechanisms (slow 
accretion, gravitational interaction, etc.) that can justify the ob- 
served asymmetric properties and are not probed by the current 
surveys. It could very well be that major mergers - that are those 
that we are constraining with our analysis - are not needed to 
explain the observed lopsidedness frequency. These surely are 
issues to investigate with future (gas multi-phase) surveys. 



5. Summary 

We have investigated the relationship between galaxy stellar age 
and stellar mass as a function of redshift in a mass-complete 
spectroscopic sample selected from the VIMOS VLT Deep 
Survey (VVDS). Up to z ~ 1.3 we confirmed the presence of 
a relationship mass-stellar ages that parallels the one observed 
in the local Universe (see Kauffman et al. 2003b). In all redshift 
bins explored, low-mass galaxies (log(M) < 10 Mq) are dom- 
inated by a young stellar population, as witnessed by the low 
Dn4000 values. For higher mass galaxies (log(M) > 10 Mq), the 
percentage of galaxies dominated by an old stellar population is 
much higher and grows regularly with cosmic time. This pro- 
cess is more efficient the higher the galaxy stellar mass consid- 
ered. This result supports the so-called 'assembly downsizing': 
the stellar population in massive galaxies formed at earlier times 
than the one in low-mass galaxies. 

We then explored SFR evolution as a function of stellar mass. 
The percentage of quiescent galaxies, as witnessed by their low 
EW[OlI] values, increases when moving to lower redshifts and 
higher masses. This trend is clearly visible both for galaxies with 
a young stellar population and for galaxies with an old stellar 
population. The emerging picture is one where low-mass galax- 
ies (log(M) < 10 Mq) are subject to bursts of star formation ac- 
tivity up to recent times, while the dominant mode of star for- 
mation for massive galaxies (log(M) > 10 Mq) is a smooth one, 
with no significant contributions from secondary strong bursts. 
Therefore not only does the bulk of the stellar population in mas- 
sive galaxies form at earlier times but the location where the 
mass also gets assembled more efficiently at each redshift moves 
to lower mass galaxies with cosmic time. 

We have also shown a mild total evolution with redshift of 
mass assembly history, in particular for the very massive galax- 
ies (log(M) > 1 1 .4 Mq) which are mainly early-type objects, up 
to the highest sample redshift (z ~ 1.3). Since the relative frac- 
tion of massive (log(M) > 10.5 Mq) early- and late-type galaxies 
evolves with cosmic time, we are witnessing a spectral transfor- 
mation of late-type systems into old massive galaxies at lower 
redshift. 

Finally, when we consider the joint distribution of stellar 
mass and star formation activity to quantify the efficiency with 
which galaxies assemble their stellar mass, we obtain a scenario 
where it is possible to account for the number of passively evolv- 
ing galaxies up to z ~ 1 without invoking any (dry) merging 
mechanism. Our observations agree well with a scenario where 
the stellar population in massive galaxies formed earlier than the 
one in low-mass galaxies, and massive galaxies themselves were 
the first to be assembled. 
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